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SU14MAEY 


With water as driving medium and delivered medium in a device sim- 
ilar to a simple Jet apparatus, the pressure and velocity fields of the 
mixing zone were explored with a pitot har; the ratio of delivered to 
driving volume ranged "between the valties 0, 1 , 2 , and 4. 

An attempt was also made to analyze the mixing flow mathematically 
hy integration of the equation of motion, with the aid of conventional 
formulas for the turbulent shearing stress, hut this succeeded on3.y ap- 
proximately for the very simplified case that a drivirg Jet is intro- 
duced in an milimlted parallel flow, while the pressure over the whole 
mix,ing field is assumed to he constant (pt, 2), 

In spite of these dissimilar assiuaptions for the theory and the ex- 
periment, the form of the measured and the computed velocity profiles 
indicates a very high degree of approximation (fig. 22). 

The pressure rise, which -vras approximated hy Flugel’s formulas, 
disclosed good agreement with the measured values (fig. 25). 


INTRODUCTION 


Surveys of turbulent mixing fields published so far deal only with 
fields of practically constant pressure (although Fortlimann (reference 
1) presents an experimental arrangement where pressure and velocity are 
measured in a completely closed, abruptly diverging channel, but fails 
to give any details about the results of the pressure measurement). 

The principal object of the present investigation therefore is to 


; ^"Untersuchung turbulenter Mischvorgange ." Forschung axif dem 
Gebiete des Ingenie-urwesens, v'bl. 12, no. 1, Jan. /Feb. 194l> pp. 16-30. 
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siirvey mixing fields with variable presstire, similar to those prevailing 
for a jet apparatus, 

A further aim was the prediction of such mixing fields on the basis 
of Prandtl’s formula for turbulent shearing stress, but tha.s was suc- 
cessful only uinder considerable limitations on accovmt of the matheTaati- 
cal difficulties involved. 

Lastly, a simplified method of calculation for turbulent mixing 
processes ira.s checked for agreement mth reality against the measure- 
ments . 


TiSEFSEIMEKTAL tElEOI) 
Exx^erimental Setup 


The size of the ax’rangement was chosen with a. view of aialring the 
survey apparatus which was used for measuring the- iiressure and velocity 
of the water serviixg as test fluid, smodl enough so that its presence 
caused no appreciable disturbance in the s.ttitude of the flow, while a 
certain minimum size- of the sux^^'ey apparatus is necessary for x-easons of 
manufacture and experimental manipulation. 

The order of magnitude of the utilized velocities is dictated by 
the demand for convenient experimental manipudation and ahsence of spe- 
cial difficulties in the positive and negative pressures anticipated 
during the measurement. Furthermore, the velocity head corresponding to 
the maximum velocity was to be .realizable within the availabDe building 
space, because it was expedient for the. a,tta.inment of a urilform state 
of flow to produce this velocity by means of overhead tank with free 
surface . 

The setup is illustrated in figures 1 and 2. A pump driven by 
electric motox* piimps the water out of the fil.led reservoir 1, through 
the suction and pressure line 2 end, . 3 into an ovex-liead tank k. An ove.r- 
flow 5 ensxires - independent of the speed fluctuations of the pump - a 
constant pressux-e head fox* the jet and for the supply of the delivex'’ed 
water, volume (secondary water), which is measured by a test nozzle 8 at 
the end of the sxipply line. The Jet leaves the driving nozzle 9 nt a 
speed of aroxuid 3.0 meters per second, mj.xes with the stirroxinding vrater 
in the cylindrical mixing chomber 13 while giving off momentum an.d 
carries it along with it. The mixture gains the outside at the end o.f 
the mixing chamber through holes 15 in the envelope of the mixing pipe 
l4, flows doTrmward throxigh a branch 16 , passes through throttle valve 17.> 
and back into resei'’voir 1. The water level in the reservoix’ is so high 
that valve IT is under water; this px-events disagreeable splashing of 
the outflowing water jet. Fox- the same purpose the level in the supp.ly 
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tank 22 for the secondary vater is also maintained high enough hy appro- 
priate' '.throttl^-Tfaive : setting to keep- the- mouth'-of nozzle -under - 
■v?ater, ,d?he kinetic energy of the outgoing secondary yater stream is 
dissipated in a small svash chamher 19. Straightener and screens 20 , 
then provide for a imiform. inflow to the driving Jet. Arranged directly 
aroimd the Jet, a n^row-cell straightener 21 (fig. 2.) with radial rihs 
ensures the elimination of circum|’erential components. To avoid rota- 
tion in the driving 'Jet the supply line 6 to the driving nozzle is like- 
wise fitted .'vrith a straightener 10. : different secondary water volumes 
are obtained Td.th different test nozzles 8. The quantity of driving 
water -is the' same for all testa. 

The measuring Instrument is a long pitot tube 23 (fig. 2) ^ carrying 
the actual survey apparatus at the tapered forward end.. Aside from the. 
dynamic pressure orifice, it is fitted with lateral orifices for record- 
ing the static pressure. Figure 3 illustrates the pitot tube and its 
particular tip. Eubber tubing serves as lead to the vra.ter and mercury 
manometers. Throttles and glass vessels acting as air chambers are 
mounted in the test lines to dampen the occasional3.y violent pressure, 
fluctuations. The pitot tube is suspended by rings 2h and, small metal 
stripe 25 (fig. 2) from a strong iron beam 26 which, together with the 
suspended impact tube, can be moved vertically up or doto by means of a 
nut and spindle 27. The height positioii for each spindle is read on a 
scale 28. The rings and strips which, on backward shifting of the pitot 
tube, came directly before or behind the test orifices, are removed. 
Distance pieces and welded-on stops 29 serve for adjusting scale end 
pointer, as well as control. The rear end of the mixing pipe is closed 
by a plate 32 through which the pitot tube is introduced by means of a 
st\iffing box, which can be raised or lowered with the pitot tube by 
hand- wheel and lead screw 33. 

Thus,; by raising and lo-wei-ing the pitot tube and by shifting it 
along its axis,, the test orifices can be placed at any desired point of , 
the vertical median plane of mixing chaiuber , with exception of the area 
adjacent to the wall, which, however, is of lees interest in the study 
of. the mixing processes. A number. of test orifices in the horizontal 
median plane along mixing pipe Ik enables pressure measurements directly 
at the wall (fig. 4). The position of the pitot tube with respect to a 
displacement along its a:cis can also be read by means of scale and 
pointer 3^. Tlie pointer is solidly connected to plate 32; the scale is 
scratched on the pitot tube. 

Calibrations 

(a) Overhead tank ..^ To calibrate the nozzle it was first necessary 
to. gage the overhead tank. To this end it -vra-s filled and then emptied 
in stages, the drained water being weighed wd the respective 
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abatements of the level measured. The mathematically obtained crossv : 
sectional area Fq of the, tank showed itsehf fairly constant oyef - the 

heighten question, .hence could be used as constant in the volume deter- 
minationt , ■ - ' ■ ' ' 


(b). Nozzles .- All nozzles are calibrated by discharge, test at 'the 
place of utilization. The necessary formulas are briefly outlined (figi 
5). (Subscript 1 refers to start of outflow test, subscript 2 to end of 
outflow test; differences in height level between start and end of test 
are indicated with . A. ) ’ 


The volume of the outflow at speed C in time interval 
equal to the decrease in volume ixi the overhead tank: 


dt 


is 


aFcdt = 


( 1 ) 


where a is the outflow coefficient. From the energy equation, applied 
toeross sections II and III, the outflow velocity- C follows at 





( 2 ): 


•where ,e is the cross-section ratio F/F^ and ^ the loss coefficient 
of the nozzle . , 


On the premise of .constant pipe friction coefficient the, energy 
equation, applied to the sections ' I and II, gives 


' ^ _ 'h 
,dH. ” ff 



(3) 


Integration 'of equation (l) by zaeans of equations (2) and (3) gives 


— ^/^ 7 ' I - 

aF 


2 ' 

e 



-J t 


N • 


I 
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Quantities t, hi. Hi, and Eg are determined hy test. Herewith the 
throughflow quantity 1^60031163 



otcFr = 




Table 1 gives the data for the throughflow quantities in the prin- 
cipal tests for the driving nozzle and the three nozzles enqiloyed for 
the secondary water nozzles-; h is the effective pressure head ineasured 
in the main test, the ratio of secondary water quantity to driving 
water quantity, which approximates the values 0, 1, 2, and k for 

the four principal test series. 


TABLE 1.- EECOEDS OF PRESSURE HEAD h, VOLUME OF SECONDARY AND 


DRIVIMJ WATEE, AND THEIE PROPOETIONALITY FACTOE AS 


OBTAINED IN PRINCIPAL TESTS 


h 

G 


(m) 

(kg/s) 


5.25 

7.19 

0 

5.28 

7.68 

1.07 

5.19 

14.82 

2.06 

4.89 

28.30 

3.94 


(c) Pitot tube .- The calibration of the pitot tube -vras effected in 
the free water Jet produced by the driving Jet, with mixing pipe removed, 
the pitot tube being suppoirted so as to swiiag in the axis of the Jet - In 
the longitudinal plane - or at any desired angle of setting with respect 
to the Jet axis. It is rotatable about its axis in every position. 

The pitot tube should have the following characteristics with re- 
spect to the recording of the two test quantities: 

1. The total head should be indicated without falsification if pos- 
sible, even in yawed flow (i.e., when the flow direction is other than 
parallel to the pitot axis) . 

2. The static pressure indication should give a necessary correc- 
tion in the simplest possible form, and be, as far as possible, inde- 
pendent of the yawed flow. 
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To meet the second z’eguirement;, it -was necessary to study several . : 
forms "before finding a suitable one. 

The first attempt was to esta"blish the limits within which the as- 
sumption, that the px-essure of the sux;rounding air prevails in the Jet, 
is Justified. Figure 6 represents two test series giving the pressure 
on the Jet axis. It is seen that the pressure "balance .is reached at a 
distance of a"bout 1-nozzle diameter (30 from the orifice. From this 
pressure it is assumed that it is equal to the pressure of the surround- 
ing air. 

All the cali"brations enumerated hereinafter were planed so that the 
momentary test orifice of the pitot tube was placed at a point of the ■ 
Jet axis 50 laillimeters distant from the mouth Of the nozzle. The meas- 
ured dynamic pressure coxrresponded very accurately to the speed computed 
for the total head from equation (2), where t = 0 and = if 62 '^ ^ 0 

for axis adjacent points. Yawed flow up to 12° and turnijag of the suxv^ 
apparatus about its axis effect no change. Thus the first of the re- 
quirements is fulfilled. 

Figux’e 7 shows the indicated static pres sxire head hg-^ plotted 
against the velocity head h^ = c'/Sg for different ang3.es p betwen 
Jet and suarvey apparatus axis (different yawed f].ows). Two cases must 
be differentiated. 

1. The flow may be diverted against the iv-ide side (a = 90°), or 

2. Against the small side of the apparatus (a = 0°) . 

The case a = 0° is obviously unfavorable for a simple correction 
formula because of the substantial depai’ture from proportionality exist- 
ing between indicated static pressure head and velocity head even at 
P = 3°. At a = 90 °, on the other hand, hgt remains unchanged 
as function of h^ as. far as p = 6°, and even for p = 9° the cali- 
bration cuxve differs but little from that for p = 6° in the rarige of 
low velocities. As assumed for the present and also. proved later on, ■ 
this Insensitivity to yawed flow extends over the greater part of the 
mixing field. In conseqxxence, the calibration ctirve for p = 6° serves 
as a basis for the interpretation. .Its equation is according to the 
measurements; 


0.007 h^ ~ hg.^;; - hg.^.^j. , 

where hg-jpijy is the true static pressxire. The other re3.ation reads 

hg = h^ + hg-t-w 
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vhere hg denotes the total head, recox’d.ed. at the tip of the pitot tube. 
Then the desired.^gLuahtiti'es h^,' and. the 'velocity c follow at: 

. : hst - 0-007 hg 

hstv = ( 

‘ ./ -' 0 . 993 - 


hd. ~ hg - hgt'Vjr (5) 


c = B hd ■ (6) 


In the course of the subsequent measurements it was found that the 
survey apparatus was apt to deviate from its calibration curve as a re- 
sult of a slight variation in surface condition (material: brass) . Hence, 
it was necessaiy to rub it slightly down with prepared chalk every few 
days and then rechalk the original curve bj'- a second calibration meas- 
urement} the process was repeated, if necessary. 


Mixing Field Sxxrvey 


In conformity with the quiantity ratios li = 0, 1, 2, and k of 

the secondary to the driving water volume, four series of tests were car- 
ried out, with about IIGO test points distribiited in the mixing chamber 
for each series. 

The survey started with the static pressure at successive points on 
the axis and on axially parallel straight lines, axially distant 20 mil- 
limeters corresponding to the identical distance from the tip of the 
pitot tube to the test orifices. Thus each position of the pitot tube 
yielded the Indicated values for total head and static pressure for two 
different points in the mlxli^g field. This fact was allowed for accord- 
ingly in the interpretation of the measurements. 

These tests on axially parallel straight lines vievQ usually placed 
at 10-millimeter radial distance partly above, partly below the pipe 
axis so as to spot any lack of symmetries. It was further found advis- 
able to deteimine the pressure and velocity distribution more accurately 
by closely spaced test points in several cross sections, particularly 
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in the part of the .mixing field adjacent, to the toivlng nozzle, where . ^ 
Telocity and pressure vaiy suhstantlalJ-y with the distance from the pipe, 
axis. 


Prior to each measurement the two. ducts in the pitot tube and the 
manometer tubing were blown out' to expel air bubbles. 

In spite of the damping deTices installed in the test lines, the 
indicated pressures fluctuated very severely. Especially disagreeable 
were the irregular, protracted fluctuations which entailed observation 
periods of an average of 3 minutes. 

Therefore the graphical representation of the explored mixing 
fields reproduced in the following should always be regarded with the 
limitation that the quantities represented in relation to the location 
are merely averages over comparatively long-time intervale. 


Interpretation of Measurements and Test Data 

The representation of the pressure field is predicated upon the 
assumption .of some chosen reference pressure, as for instance, the 
pressure ,,'5fhi. in the annular slot at the driving nozzle entrance (fig. 

2) .. With’ h'g, and ..h'g^. denoting the read- off ’;ra,ter column, levels of . 

the principal test corresponding to the pressure at the tip and at the' 
orifices of the survey apparatus (all levels h being measured with 
respect to atmospheric pressure), there is obtained 


hat h'gt - hi; hg - h'g - hi 

Pressure and velocity are computed by equations ( 4 ), ( 5 )^ and (6) and ' 
then plotted against the distance from the driving nozzle entrance for 
equal distances from the axis. Then a number of axially normal sections 
are traced through the' mixing chamber and for these also the measured 
pressure and velocity distribution is plotted point by point, s,pocial 
attention, being paid to the cross' sections with the closely spaced test 
points., , " ' ... 

In the following the pressure corresponding to the value hstw is 
is indicated with p (p = 7 hgtw) • Figure 8 shows the pressure and 
velocity • in two cross ■ sections, .chosen at random. ■ 

The next problem is to. complement the . pressure and velocity dis-: 
tribution established, point by point 1)7 a .continuous .area.. This makes ■ 
it necessary, to fix a rule by which certain irregularities of the 
measured curves are smoothed. 
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At greater Tallies for ratio p. it was found that the results of a 
number of the surreyed cross sections were not exactly symmetrical with 
tho ’axls Of the'lnlxliig pipe I The ratios -vreire plot^^ at exaggerated 
scale, as in figure 9 , where TDC is the geometric axis of the mixing 
pipe, with the points Oi, Og, O 3 , and so forth; The points 1, 2, 3t 

1 ’, 2 ' , and 3 * are intersection points of the surveyed velocity curves 
with parallels to the axi s of a bscis sa. The points Mi, Ms, and M 3 

divide the distances 11', 22';i aiid 33' ' in eijual parts . . For the fur- 

ther Interpretation the line through Mi, Ms, and M 3 , and so forth, 
was hent straight and regarded as axis. The asymmetry is, moreover, 
slight j the distance of this axis from the line AA is, at the most, 

1.5 millimeters. The cause of this irregularity may lie in the asym- 
metry of the inflow for the secondary water; ^nor defects In the shape 
of the mixing pipe may contrihute; a deflection of the pitot tube would 
also act in this direction. 

The asymmetries existing in the pressiAre curves were smoothed in 
the same manner. In the area adjacent to the nozzle the measiired pres- 
sure dlstrihutions over cross sections frequently showed a form as re- 
produced in figui« 16 . The measured minimum pressures are not aliSte hy 
otherwise pressure curve symmetry . In. such cases the arithmetic mean 
value serves as presswe peak for the interpretation. The cause of this 
iri^gtilarity inay, aside from irregular inflow, he due to the scattering 
Of the test points, which is comparatively great at these very points. 
The explanation for this unusual pressrire mliiimum will he given else- 
whero . y- ' ‘ 

' Figures 11 and; 12 represent sections thfqu^ the pressure and 
velocity areas for mixing ratio p. » 2 ; figure 13 , the preSsuire die trl- 
hution in individual sections, aged figure l4,the velocity profiles for 
p » 2 (x « distance from driving nozzle) . Once the longitudinal and 
cross sections are definitely fixed> these two Sets of curves define the 
pressure and velocity areas. For lucid representation these areas are 
Intersected hy a series of parallel planes, the analytical expression of 
■vdiich is p * constant, and c * constant. The section curves eire 
represented in figures 15 to I 8 .; The streainli^ deteimined hy 

graphical Integration on the basis of the measured velocity profiles. 

The streamlines are dOteimLned on the simplified assumption that 
the measured velocity, with sufficleht accuracy, shall he equal to its - 
axial component (c K ^ax^ * points of maximum inclination of 

the streamlines the difference between velocity c and axial component 
Car is- 3 *0. k percent^for p: » Oj-'considerahiy less for p * 1, 2, 4. 

The secured streamlines are also shown in figures I 5 to 18 . 

Error in the Static PressTire Record at Points of Great Velocity Gradient 

Both Schlichting (reference 2) and FOrthmann (reference 1) have 
already pointed out that the static pressure readings obtained in flows 
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with prohoiinced lateral velocity field .are defective.: >13ie principal 
cause is most like Ij" to. he 'fotmd in the 'Strong -ptxlsating transverse -vj: 
velocities: accompanying such flows * as is' evident froia the •■Galihration:':.,- 
curves in figure 7, the survey apparatus is not, .insensitive to this yawed 
flow. The static pressure record is smaller than assumed in'the, utili^^ 
calibration curve. •’But it is : somewhat amusual that, according to’ ', 
Forthmann (reference l) , these negative pressures can amount to as •.much;, 
as 15 percent of the impact pressure of the' maKimum speed, in the test'.'..;, 
cross section,^ as against only ahoul;. 1 i>ercent in the present case . , 
Approximately the same value was obtained with a survey, apparatus) of , ' ' 
circular cross section 6 millimeters, in diameter arid lateral test ofi- . 
fices., • „ 

. Figure 10 shows the "form of tlie' discussed 'negative pressure maxi- 
mums. it most likely is a defective raeasureme^i-fc , and for •the following . 
reason; since the streamlines in this area are nearly parallel and 
straight, obviously,, the pressure on normals can vary but lit-tle . . The 
pressure measxxrement on the axis is correct or practically so, because 
the interchange of masses is very small. The probable pressure distri- 
bution is therefore somewhat as indicated by 'the dashed curve in- figure 
10. . This probable pressure distributlon.in the septate, sections Is. 
plotted so as to form a. continuous surface (dashed curve in figs., 15 to 

Experiments with Driving Nozzle Extending Deeper into the Mixing Chamber 

The position of the driving nozzle relative, to the inlet of .the 
mixing pipe may raise the suspicion that wi'th the chosen arrangement, the 
total rise is adversely affected, insofar as tb®. .S'^art of interminglii^g. / 
occurs in an area where certainly , some irregularity of flow.still ex- 
ists. It therefore seemed .desirable to' measure ' the wall pressures with, 
a nozzle extending farther, into. the miming chamber in-order to detect , 
any .eventual variation of •the pressure .field. , ' • 

Figure 19 shows' the .lengthened and the normal, nozzle.. F^jr.the;"/" ,, 

^It is not gtii-be clear from the cited reference whether idiis 15 
percent refers tp the local impact pressiure or to that of the .maximum 
speed in the test cross section. The second interpretation is rather 
the logical one. ^ . - 

^The greatest negative' pressure peaks , occur in cfpss . sections , in 
which the core. of the, driyihg Jet has not ah yet ' di'sintegratedj the max-, 
imum speed in these sectione’ is therefore idehticaJL ■with the discharge 
velocity of the driving Jet from the nozzle. 
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tests with the normal nozzle the pressure 7hi in the annular slot of 
the nozzle (fig. 2a) aerTed as ref erence j)^^ssure . B^e^ pressure 

ahd'pressure 'Thg (test station, fig. l) tfie following relation exists 

7(ha - hi) = a 


a assuming a specific value for each mixing ratio. Now, however, this 
pressure cannot he measured directly and is therefore calculated from 
the above relation; a is taken from the previous tests, hg is read 
off, as usual. 

The tests indicate that the presfsiire rise with long driving nozzle 
is very similar to that obtained with short nozzle , except for being 
shifted backward for th© excess length of 120 millimeters of the long 
nozzle. The somewhat unfavorable inflow for the start of the mixing 
with short nozzle has therefore no adverse consequences on the total 
pressure rise . 


THEORY 

Simplified Assumption and Method of Solution 

An attempt was made, at first, to make mixing processes of the 
kind seizing as basis of the experimental study, amenable to theoretical 
treatment by means of Prandtl’s shearing stress formula, but these at- 
tempts were frustrated by mathematical difficulties . Then it was 
attempted to make flows which approached the explored mixing flows as 
closely as possible, amenable to theoretical study by simplifying as- 
sumptions. 

In the following the intermingling of a Jet with a parallel flow 
acting in the same dinrection by axially symmetrical arrangement and the 
assumption of constant pressure over the mixing field is treated as a 
flow of this kind. As for the previously computed mixing flows (refer- 
onces 3, 4, 2, 1'), the follpwing assumptions are mde: 

1. The mixing processes in cross sections with varying distances 

from the start of mixing shall be geometrically and mechan- 
ically similar. 

, 2 . Ihe mixing path Z shall be constant over a cross section and, 

therefore, proportional to the. width of the mixing zone^ fer 


^Assuurotions 1 and 2 were originally made by Prandtl (reference 


5 ). 


12 


NACA TM No. IO 56 

different cross sections, according to the first assumption. 
The procedu:pe is hriefly as follows: asstimptions 1 and 2 in 

conjunction vd-th the impulse theorem themselves make it pqs- 
sihle to find a 3.aw fbr the increase in width of the mixing 
field (as in references k and 2 hut originally shown in 
reference 5 ) . Herewith is suggested a formula for the ve- 
locity, which, together with Prandtl*s shearing stress 
formula, reduces the employed simplified equation of motion 
to a common differential equation, the soliition of which 
presents no difficu3.ty, 

So3.ution 

(a) The Prahdtl shearing stress formula .- In side-'by-side flow of 
layers of density p, with a velocity, gradient largely perpendicu3.ar to 
the principal velocity, the interchange of momentum introduces an appar- 
ent turhulent shearing stress, which is expressed hy 


T 



Sy 


) 


du’ 

Sy 


( 7 ) 


where u’ is the average local velocity in x-direction, y the coordi- 
nate at right angles to it, and I the convection path, The average ■ 
local transverse velocity v' of the mixing process, which causes the 
momentijm interchange, is proportional to 


hence 


Su' 




(8) 


(h) Formula for the velocity. - From the similarity of the state of 
flow in different cross sections follows 


I ~ db 


( 9 ) 
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and for slinllarly placed points of different cross sections the' propor- 
tionality, relation 

. . •s*' - 


a is a proportionality factor, h the vidth of the mixing field, u 
the relative velocity in respect to. the parallel flow of velocity U; 
hence 


u‘ = U + u and whence 


Su’ _ du ’ 

Sy , ' Sy V , 



With it 8 oid equations ( 8 ) and ( 9 ): v'-^ou. In addition 


, flh 8h Sx 

dt dx dt 


u' 


8b 

8x 


(U + u) 


■ 8b 
8x 


The only case considered hereinafter is that where u is negligibly 
small relative to TJ, which is at lea,st admissible for sufficient dis- 
tance from where , the mixing starts. Then 


8b 

v' = U T— ~ ecu 
8x 


,db. 

or IT — r odx 

' , ' ' • U'' 


( 10 ) 


A. relation between u and b'"'-fbr integrating equation (10) is afforded 
by the impulse formula, applied to the control space according to figure 
80: _ 

b/-s ' " . 

-i— = 2U / uydy 

2lip.^_ J- 

o 


( 11 ) 


wherein one term u^ is neglected against 2'Uuj J is the in^julse of 
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the driving Jet,, only the increased velocity relative to the; parallel 
flow coming in consideration as velocity portion. The integra].- is pro- 
portional to a product uh^, when u is the value for different cross 
sections taJcen at similarly placed points. With it xx'-xi/h^ {m^/sec) 
n "being a constant. The integration of eq.uation (lO) then gives 


■ ' -3^ " 

X ~h - or 

n 


where s = (n/U)^/^. for short. Similarly placed points of different 
cross sections therefore are situated on curves of the family 



y y 

Tj =: _ j= — ~^-=: constant (12) 

Tu gj.1/3 


The velocity formula suggested in this and sim3.1ar cases (references 3^ 


• ' u = Usq)(x)f('n) - ( 13 ) 

where q) is a function of x and f a, function of gives, in con-. 
Junction with the condition that according to ectuation (li) the inroulse 
in ..each cross section has the same value, hence must he independent of 
x: 


cp(x) x^/^ = constant 


Herewith eq^uation (13) is transformed as follows (the constant being now 
included in f(q)): 


- 

. •y-2/3. 


f(ll) 




U 


(14) 
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(c) Integration of equation of motion .- .Figure 21 represents a 
rii^ element of the explored flow which, is hounded hy stream surfaces 
and planes noraal to the x- axis. The impulse formula for this element 
under the assumption of constant pressure reads: 


dR - dB = 0 


(15) 


dR signifying the increase of the frictional or shearing force in radi- 
al direction, dB the increase of piomentum in axial direction. Dis- 
counting the minor slope of the streamlines relative to the axis, both 
expressions are developed by the continuity equation and the relations 


dB =-^ dx; 
5x 


dR = — dy 

oy 


where B = me (m = specific masa of flow, c - its velocity), 

R = TRitydx. So, finally, the equation of motion in "natural” coordi- 
nates reads 


5c 1 5(jy) 

5x py 5y 


(16) 


In view of the nearly axially parallel flow, c ~ u' can be reasonably 
approximated to c » u' — U + u. This affords, ■VTith the effected omis- 
sion of u relative to U: 


^ _ 1 5(Ty) 

5x py 5y 


(IT) 


The expressions for u and T , according to equations (7^ a^d (l^), 
are now written in this equation, whence, with the use of equation (12) 
and the relation ; . : 



or 


I = ksx^/® 


(18) , 


I 
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where . k denotes an as yet unknown quantity. 


— [2f(n) + f (ti)ti 
3 Sq 


or 


(19) 


(2r,f + Ti2f’)dri = 3 k2d(Tjf'2) 


f "being the derivation of f with respect to t). The eqmtion can he 
integrated once; 

Tl^f _ 3k2T^f'2 „ c (20) 

For Ti = 0 follows C = 0, since f and f remain finite; 

T]f - 3k’^f'2 = 0 (21) 


Equation (21) can he integrated after separating the variables; 


f = 



( 22 ) 


The constant Ci follows from the houndaiy condition that for u = 0, 
q itself has a well defined value For u=0, f(q) is obviously 

also equal to 0, hence 



3^3^ 


(23) 
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Eqtiatlon (22) then reads: 


f(il) 






(24) 


A further houridary condition defines the still unknown qtiantity s nnd 
the liinitii5g value t]q: the difference of the impulses over each ay.ial3y 

normal plane upstream and downstream from the: driyipg Jet entrance must. , 
yield the impulse J of the dfiVin^‘ according to 'Sq,iiationB 
( 11 ), (l4), and (24): 


J 







dTi 


(25) 


The impulse J and the , temporarily ;cgadeterminod s are constants. 
Since tj and k are dimehBion3.esB the formula can he written 


J = pU“s3 


This necessarily^ specifies compliance vrith the condition 


1 = 


4rt 

oyjk)' 


ho 




s/e 


TJ 


0 


■whence : 


= 1.756 k^^. ' 

and, •with 'the use of equation (23): 

Cj. = -0,448 


^ ^ (25a) 
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Then the complete exparession of the velocity u is: . 



The hovindary of the affected zone has, according to equation (25a), the 
form 


y^ = 1.756 k (27) 

^ pU2 ^ 


From equation (26) the velocity distribution in a section perpendicular 
to the axis then follows suitably dimensionless as quotient' u/ug,, Ug 
being the velocity in the axis. Put for the sake of simplification, 



3^3/3 


the velocity on the axis is 


ua = HCi2 

ft-ntS put, furthermore, tj = St)q, where 5 may assume any numerical value 
from 0 to 1, there is obtained 

L = (83/2 - 1)2 
% 
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The computed velocity profile is represented in figure 22. 


COMPASISON OF; CALCULATION AlfD EXHSEIMENT 
Decrease of Velocity on the Axis 
and, Increase of Width of the Mixing Field 

•According to equation (26) •the velocity on the ctirves ■ tj = constant 
decreases donfosmiahly to a power law. This decrease should he particu.- 
larly easy to verify for the axis (ri = 0) . According.to a similar, 
simple power law the width of the ioixing zone increases with the di. stance 
from the point where the mixing starts, equation (27). Both laws are 
developed on the premise of constant pressui-e in the mixing chamber 
without considering boundary walls . 

Because of the substantial pressure rise occurring in the measured 
mixing field, the veloci^by must naturally drop much faster than obtains 
by a calculation on the assumption of constaait pressure, and the in* 
crease in mixing zone -width in the cylindrical pipe also will differ 
from -that in the unlimited parallel flow. Therefore it is not. a. priori 
possible to find these theoretical laws realized in the executed mea.s- 
tirements. ■ And in consequence it is also impossible to define from these 
tests the theoretically inde-terminable . cons-fcant 'k empirically in its 
order of magnitude. 

For the case of mixing of a round Jet with the surrounding stil3. 
air the ratio, can be taken from the tests wi.th the value O.O738, 

according to Tolimien (reference 3) • On the assumption that I/Jq 
shall be 0.0738, equations (12), (18), and (25a) give the value 
k 5:1/30. 


Form of Velocity Profiles 

Although the assumptions for theory and test differ, the form of 
the measured velocity profiles suggests a compad*ison with -the theoret- 
ically computed velocity distribution. Of cdxirse, this involves only 
cross sections of the mlxiiig chamber in which" .. 

,,, 1,. The,,so-imd, pf the driving Jet is already com- 

pletely gone and for which 

2. The mixing motion has not advanced too close to the -imll 
proximity, so that it is still logical to define a 
theoretical inflow velocity U of the delivered water: 
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U =- 


r(Fr - ^’o) 


G -2 (kg/sec) is the delivered q.i^Jitity of water, the mixing pipe 

cross section, and F^ the driving Jet cross section at the exit fromi 

the driving nozzle.. Figure 22 shows, aside from the cou^juted velocity 
distribution, several of the measured velocity profiles for n = 1,. It 
is seen that an acceptable agreement exists between theory and test 
throughout the particular zone, in spite of the fact that the theoreti- 
cal assumption that should be small compared to U, does not hold, 
for the cross sections employed for the coiapariBon. 


APEROXIMATE ANALYSIS OF THE MIXING PROCESSES 
Fliigel’s Method 


In Fii^el's article (reference 6) on the calculation of Jet appa- 
ratus, a dragging force effective between driving Jet and entrained Jet 
is presented in an unusually simple manner. Figure 23 shows a driving 
Jet which entrains the surrounding fluid through mixing motions. For 
simplicity it is assumed that in each cross section, the velocity Cg of 
the delivered medium and the velocity cx of the driving Jet shall bo 
constant. The effective dragging force dE on a surface portion dO 
of the boundary area between driving and entrained Jet .shall then be: 


dE (ci Ca)^ dO 

■ ,2g 


X is a kind of friction coefficient. Next, the Impulse formula is ap- 
plied to a small sector of the driving Jet and the delivered Jet. After 
various transfonnations a differential equation is obtained which, among 
others, is solvable also for the assumption of constant throughflow 
cross section F and thus makes it possible to compute the velocity 
ratio Cs/ci; that is, the degree of mixing for different distances froa 
the start of mixing. The final formulas according to Flugel's report, 
are briefly introduced: 
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G-sCix 

CeX = — 

F/Cix - Gi 


(29) 


G"3_ 

Ap = — (cio - Cix) 
gF 



Cso) 



Geo Gsx ^ 


(30) 



Herein: 




°130 

driving Jet velocity 

1 

at 

reference cross section A-A 

G20 

velocity in entrained Jet 

j 


(fig. 24) 

GlX 

driving Jet velocity 

1 

at 

distance x from the reference 

®2X 

velocity of entrained Jet 

J 


cross section 

Ki 






Ks, = ysTT/ / I (_A!._ ) 

V Gi + G2 X - pjjj^ 
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Pm 


~ J PoPx 

Ct 2 + Crj_ 

Frcio 

Gi + G-2 






£L]_ = kx (l + ^2) 


ki 


3»2> ^ 


(ki - ka) 


a, = 


= 


ka 

"i" 

k 


ki = 




^■1 + ^2 
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G-2 

Icg = 

Gx + &2 
kg = kx k 2 

Gx quantity of driving water 

G2 quantity of delivered water 

Ap pressure rise relative to pressiire in reference cross section 
F cross section of pipe 
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X . coefflcierLt. of f fiction, for the pipe ’mi.ll (assiJiaod for the nxansr- 
ical calculation at X = 0.005) 


Verification of the Dragging Jorce of the Driving Jet 

Tho survey 'Of : the velocity field and the- pressure field in the mix- 
ing chamber can serve as- -basis for an empirical prsdiction of the, 
friction coefficient from equation (28) . Hio calculation is espe- 
cially simple for, mixture ratio p,= 0^ the quantities and a4 of, 

equation (28) are equal to zero in this case. Tho velocity 

: cix-= C 2 X = G1/7P 

in .the ■ smoothed cross section is then computed and from the representa- 
tion of the- pressure and velocity field the distance x from the start 
of :mixing where velocity and : pros sure may be regarded as smoothed, - is 
•ascertained. Every thirg in equation (28) is then kno-wn except the - • - 
looked for x • • . . ' 

•' A similar attempt for tho case of n =* 1, 2, and 4, encounters the. 
difficulty that tho expression Pq Cio/cix the smoothed cross section 
becomes equal to unity. Then the third and foxn-th term of equation (28)' 
become infinitely great; that is, the equation states that tmiform ve- 
locity occurs only foi- infinite pipe length. 

To determine x in. these cases, select mixing field sections of 
not quite completed velocity compensation and read off an average veloc^ 
Ity value for the area of the driving Jet from the surveyed velocity 
profile, which then, entered as c^x^ together with the distaiace x of - 
tho particular cross section in equatioii (28), makes it possible to com- 
pute the value x> choice of cross sections should not go too far in 

the area of the almost smoothed zone, else the expression po Cio/°ix 
approaches unity tOo much and the denominators of the' third and fourth 
ten® become almost zero; hence an inaccuracy of the read-off term acts 
comparatively very falsifying bn the result . 


Tabic 2 contains the computed x i^ reference sections,' each' 
for the foxir- mixture ratios. 

: TABIE’’ 2. - VALDES' x ..*^0 BEFTKEITC^ '' ' 


l-l 

0 

1' 

2 

4 

X (m) 

::o: 95 ' 

■'^.;8o;' ^ 

. . - , 

" 6.90' 

1.60 ' 

'1.15 ' 1.35 

cix (m/sec) 

.30 

■ T^lO • 

• .80 

■■' 1.40 

1.15 ■ 

1.70 1.50 

X 

.102 

.093 

.119 

.087 

.108 

.10 ' .127 
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So for the use of equation (28) a value of x ** shotild he reason- 
able . • 


The Pressure Pise . 

The pressure rise is preferably obtained by calculating the respec- 
tive value X from equation (28) for a number of the assumed values 
Cio/cirJ the value x =0,1 is utilized. After 'O^x is computed by 
equation (29) the respective Ap is then found-from equation (30). 

Figure 25 shows the pressure distribution measured along the wall 
of the mixing pipe and that computed by the cited method. As regards 
the total pressure rise, good agreement obtains between theox-y end test, 
expept for the fact that the measured pressure curves are considerably 
-shifted into the mixing, chamber as compared to the calculated curves. 

For the case p = 0 the theoretical pleasure rise is slightly less than 
actually obtained by the measurement. This might be due to the fact that 
for this case an exceptionally heavy vortex ring is foirmed in the for- 
ward part of the mixing chamber, the forward £ind backward flowing part 
of which contributes impulses in the sense of an increased pressure rise 
.(fig.. 2 k) . The same conditions are to be found to a lesser exteixt for 

The pressure rise computed for p, = 0 presents another unusual 
feature; that is, the curve Ap = f(x) (see equation (30)) exhibits a 
vertical targent at point x = 0.978, or, in other words, the cited 
x-value represents a regular maximum of the curve x = f(cio/cxx) (see 
equation (28)) at the point Cix = G/(7F). This part of the pressure 
curve is reproduced at the right-hand side of, figure 25. The dashed 
branch of the curve, for Cxx < Gr/yJ has no longer any pijysical meaning. 

The pressure measured along the mixiig pipe wall represented in 
figure 25 is not, as heretofore, refei’red to the pressure ^ yhi at the 
driving nozzle entrance, but, like the calculated pressure rise, re- 
ferred to the pressure p* in the reference section A-A (fig. 2 k), 

A-A was placed at the start of the cylindrical mixing pipe to ensux’e 
the simplest possible application of Cso- i'he driving nozzle, may be 
visualized as beiig lengthened up to the section A-A, so that for the 
inlet impulse at velocity cio, ’ the (numeri calls’- very small) increase 
due to the pressure drop (yhi - p*) must be taken into consideration.. 

The determination of p-^ proceeds from the recorded pressure yh^; 
the equation of energy gives the pressure drop 'kke veloc- 

ity increase from test point to inlet cross section. 


. Hence : 
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Ap2e = rha-l-0'.E^5#..lT-r- (c*® - Ca®) 

2g 

With the assimption that the mixing pipe section through contraction as 
a result of the somewhat unfayorahle- form of inlet, as -^rell -as the ex- 
istence of the ring vortex, is hlocked hy ahout l/3, the follomng 
values for Ap2e are obtained: 

d . 0 1 ' 2 

Apa©,.; ;• 0 8 •• 32 ; V- ' 116 f -mm of water , column 

Then, if the recorded pressure difference /(ha hx) for each mixing . 
ratio is denoted by a, the difference between tlae previous apd' the 
presently employed reference pressure is ^ 


7hi - p* = Apse - a 


The subsequent compilation contains the tfest values a and the desired 
differences 


1^ 

0 

■ r ■■ 

, , . 1 

2.‘ 

4 — 

a 

0 

. . ^ 

lA;.;' ■ 

26 mm j water 

(Apae - a) 

0 

4 

22 

? column' 
90 ram 1 

t 


Hence, to arrive at the pressure rise plotted in figure 25 j the value 
(Ajtee - a) must be added to the px-essures shown in figures 15, l6, IT, 
and l8. ' ’ ‘ 


Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 1,2, 2a, 3, 4 



i, tank; 2, intake line to pump; 3, pressure line from pump; 

4, overhead tank; 5, overflow; 6, downpipe to driving nozzle; 
7, downpipe for secondary water; 8, test nozzle for secondary 
water; 9, driving nozzle; 10, straightener ; 11, stop valve; 

18, stop valve; 13, mixing chamber; 14, mixing pipe; 15, ori- 
fices in mixing pipe (fig. 2); 16, deflection branch; 17, 
throttle valve; 18, inlet from water piping; 19, small swash 
tank; 20, straightener and screens; 21, straightener with 
radial ribs (fig. 2); 22, inflow tank for secondary water; 

23, pitot tube (fig. 2); 84., suspension rings; 25, brass 
strips; 86, iron beam; 27, nut and lead screw; 28, scale and 
pointer for reading level; 29, stops for control of level; 

30, tension wire; 31, tension weight; 32, plate with stuffing 
box; 33, hand wheel for raising and lowering plate; 34, point- 
er for axial displacement of the pitot tube. 



Figure 3,- Pitot tube and tip. Figure 4.- Disposition of test 

orifices . 



water column 
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Figs. 5,6,7,8,9,10 




fin W W 


Figure 5.- 
Calibration 
a St -up. 


a-30“_ 

0 .- 0 ° ~ 


c-i7m/stc 


c-Tm/sec 


Distance from nozzle entrance 

Figure 6.- Static pressure 
reading at various 
distances from the nozzle 
entrance . 
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Figure 7.- Static pressure reading 
plotted against dynaunic 
pressure at different yawed flows 
(P = angle of flow, oi =■ angle 
torsion about axis of survey appa- 
ratus ) . 



Figure 9.- Schematic representation 
of the asymmetry at 
‘measured velocity profiles. 


Figure 8.- Velocity 
and pressure 
profiles for jji = 1 
and p = 4. 


fXJ ^ 0 Zt 


Figure 10,- Pressure profile near 
driving nozzle; the 
probably correct pressure is ^ 
indicated by dashes. 
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Figs. 11,12,13,14 



1300mm noo 


Figure 11.- Velocity distribution over tbe length of the 
mixing field at different distances from the 
axis (section through the velocity areas) for ii = 3. 



Figure 13.- Pressure rise over the length of the mixing 
field at different distances from the axis 
(section through the pressixre areas) for |jl = 3. 
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' ' i ' ' iL Figure 14.- Velocity profiles 

^ in the individual sections of 

Figure 13.- Pressure distribution "tlie mixing chamber for \x - 2. 

in individual sections of the 
mixing chaiaber for p, = 2. 
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Figs. 15,16,17,18 
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Streamlines 


Lines of equal velocity (cax •const, 


Lines of equal pressure (p'ConstJ[kglm,'J 




'KO-7iO -m 


Figure 15.- 


Streamlines 


■too 200 300 wo . son soo too Si.. 

Lines of equal velocity (car -const) [inlsec] 


Lines of equal pressure ip -const)! kg! mt] 


Ljnes_ of equal velocity fcai-const)[m,lsec 
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Lines of egoal pressure (p^.const)[hglTn^] 
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Figure 17.- | 


Streamlines 


Lines o f eg oaJ pressure (n ‘consf)[kg!m/] 


Figures 15 to 18. 
for different la.. 


Figure 18,- p, = 4. 

Pressure and velocity distribution over the 
mixing field and variation of streamlines 


























Figure 23.- Schematic sketch of driving jet, which entrains 
the surrounding medium as a result of mixing 
motion (from Flugel report). 



Figxire 24.- Ring vortex with impulse oontrihution in the 
sense of increased pressure rise for in and 
outflowing portion. 



Figure 25.^ Measured and calculated pressure rise. 
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